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ABSTRACT 


Measurements of the sound radiation patterns in water from 
the flexural vibrations of a clamped-edge steel disk have been 
made and are compared with the results of theoretical calcula- 
trons made by Alper and Magrab (Journal of Acoustical Society 
Of America, Vol. 48, Number 3, pp. 681-691, 1970) for the two 
lowest order circularly symmetric modes of disk vibration. 

Although some differences were expected and were found due 
to the experimental condition which only approximated the assump- 
tions made in the theory, it was found that the major features 
of the measured patterns agreed reasonably well with the theo- 
retical pattern. 

The results have applicability to the design of sound 
source which could be used in underwater tracking of vehicles 


On a test range. 
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Νους τον 


A high frequency broad-band transducer which has a broad 
radiation pattern in a half-space is required for acoustic 
tracking of underwater vehicles in situations where the vehicle 
carries the sound source and the receivers are fixed. [1] The 
underwater ranges at the Naval Undersea Warfare Engineering 
Station (NUWES), Keyport, Washington, are of this type. During 
tests of production torpedoes and new torpedo designs, the tor- 
pedo is tracked in three dimensions as it moves through the 
water. 

The tracking is accomplished using an array of hydrophones 
mounted on the bottom of the acoustic range which receives 
acoustic pulses radiated periodically from a projector mounted 
on each vehicle on the range. Currently, an acoustic projector 
mounted flush with the bottom of the torpedo sends out a 75kHz 
coded pulse that is picked up by the hydrophones. A shore- 
based computer receives the information from the range hydro- 
phones and computes the location and speed of the torpedo 
Muring the run. 

For the NUWES Keyport ranges, the signal to noise ratio re- 
quirements and the expected maximum transmission loss between 
source and the receivers places a minimum requirement on 
Source level which is quite high. In addition, the pitch and 
roll motions of the torpedo and the effect of surface reflec- 


tion of the sound has led to a specification of a broad beam 


EI 





EN Specifically, current requirements are for a uniform 
source level within a conical volume about 150? wide centered 
on the vertical axis of the transducer. 

The projectors are required to operate acoustically un- 
loaded (in air) 95 percent of the time, since the unit is | 
turned on before launch and continues running after recovery. 
Another requirement placed on the projector to be used in 
higher speed vehicles is that the radiating surface be flush 
with the body of the vehicle. 

These various requirements place extreme demands on the 
physical properties of materials used in the transducer and 
create difficult problems of building a transducer which also 
has reasonable life expectancy and reasonable cost. 

Another approach to this problem is the use of a flexural 
wave transducer as a broad-beam, broad-band-width source which 
would be readily adaptable to flush-mounting, a necessary re- 
quirement for high speed torpedoes. The possibility of using 
the torpedo hull itself as the radiating surface is especially 
attractive for the low drag, laminar flow shapes expected to 
become used in the future. Present designs use a transducer 
which is inserted into a recess in the hull. Although this is 
relatively smooth, it is not adequately smooth for ultra-low 
drag laminar flow. 

Some early work in this area was done by Barlow at NUWES 
(private communication). The first previous work at the Naval 
Postgraduate School on this subject was done by LT Sevdik [1]. 
His calculations were based on classical plate theory as de- 


Scribed in Morse and Ingard [2] and Malechi [3], which does 
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not take into account of effects of transverse shear and 
rotatory inertia. Among the results of Seygdik's work was the 
noting of damping of the plate natural frequencies by fluid 
loading and of inconsistencies in the change of flexural wave 
phase speed with frequency for transducers of differing 
thickness. 

Mindlin [4] includes these effects in a three-dimensional 
plate theory and shows that classical plate theory begins to 
differ markedly from three-dimensional theory at frequencies 
above that at which the flexural wavelength is about ten times 
the plate thickness. His theory gives flexural-wave speed 
values for a given ratio of plate thickness to flexural wave 
length which are lower than those from classical plate theory. 

Feit [5] considers the case of an infinite plate excited 
by a simple harmonic point source and which is fluid loaded on 
one side. He discusses the significant changes in the theoreti- 
cal radiation patterns that result at frequencies above and 
below the critical frequency, that frequency at which the 
flexural wave speed in the plate equals the speed of sound 
in the fluid. His work is based on the Mindlin-Timoshenko 
plate theory. 

An experimental investigation was conducted of the effects 
ОҒ fluid loading on axially symmetric flexural waves ina 
ENEcular aluminum plate 25 in. in diameter and 5/16 in. thick 
by LT Jarvis [7] in 1977. Jarvis' results are in good agree- 
ment with the Mindlin-Timoshenko plate theory for the speed 
Of the flexural wave as a function of frequency. Jarvis’ re- 
sults clearly show a change in the radiation resistance as 
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frequency increases, as evidenced by the increased damping of 

the plate natural frequencies with increasing frequency. At 
frequencies below the critical frequency, the fluid adds iner- 
tia which causes a lower flexural wave speed. Above the critical 
frequency, the wave speed increasesabruptly to about the value 
for the unloaded plate, indicating that the radiation reactance 
(inertial) becomes very much smaller above the critical fre- 
quency. 

Until the work of Alper and Magrab [6], a "closed form" 
nOn of the coupled vibration of a finite circular plate 
clamped in an infinite, rigid baffle had not been obtained for 
either classical plate theory or the Mindlin-Timeshenko plate 
theory. 

The first treatment of this problem was given by Lord Ray- 
leigh [8] Mindlin [10] andDeresciewicz [9] have shown that the 
normal modes and natural frequencies of a plate of finite 
thickness are different and more numerous in a given frequency 
band that are given by classical plate theory. 

Reisman and Greene [11] have studied the "in vacuo" response 
of clamped circular plates to suddenly applied loads using the 
Mindlin-Timoshenko plate theory. They found that except for 
the first three or four natural frequencies of very thin 
plates (diameter to thickness ratio 40 or more) there was no 
Similarity between the natural frequency spectra of the Mindlin- 
Timoshenko theory and classical theory. The difference in 
transient responses were also significant. 

In the Alper-Magrab paper the acoustic pressure distribu- 
tion is obtained in a half space bounded by an infinite rigid 
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baffle containing a clamped circular plate which is driven 
from the "in vacuo" side by a simple harmonic force having an 
тату spatial distribution. A closed form solution is 
obtained which is valid for both the near and far fields. The 
analysis uses oblate spheroidal coordinates to obtain the solu- 
tion of the wave equation in the fluid. These coordinates are 
chosen because a bounding surface of one of the coordinates is 
represented entirely by the circular plate, and a bounding sur- 
face of another coordinate is represented entirely by the in- 
finite baffle with a circular opening the same size as the plate 
diameter. The third coordinate is identical to one of the 
polar coordinates, the coordinates in which the plate motion 
is obtained. As a result, the otherwise fairly complicated 
coupled problem is reduced to two simpler boundary value pro- 
blems which can be solved separately and the solutions coupled 
by the continuity condition at their common boundary. The re- 
sulting solution for the acoustic field is given by an expansion 
in spheroidal wave functions. Pressure at any point in the 
Field can be evaluated by inserting the value of the wave func- 
tions at the point. 

This is in contrast to previous known attempts to solve 
this problem, all of which have resulted in the need to perform 
difficult double (surface) integrations numerically with each 
integration valid for only one point in the field. Furthermore, 
all Known previous work has been restricted to the classical 
plate theory. 

In the following pages, a brief review is presented of 


Alper's theoretical work as reported in his doctoral thesis [16] 


185 





and from the point of view of using the results as a design 
tool tor a transducer. 

As part of the present thesis, experimental work was conducted 
as a test of the Alper-Magrab theory, using measurements of the 
radiated beam pattern from a transducer. These results are 
presented and compared with Alper's calculation. 

BS, a modification of the Alper-Magrab computer pro- 
gram is included as an appendix in a form suitable for running 


on the Navy Fleet Numerical Weather Facility Computer (CDC- 


ESUD-' HAL"). 
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fia  convpOkel [CAT CONSIDERATIONS 


A. APPROACH 

An outline of the approach to the problem [16] of calculat- 
ing the sound field from the flexual vibration in a clamped 
edge plate is: 

(1) The wave equation governing the motion in the acoustic 
fluid is transformed into oblate spheroidal coordinates. 

(2) Using the Mindlin-Timoshenko plate theory, the eigen 
frequencies and eigen functions of the plate are determined 
for a clamped circular plate using polar cylindrical coordinates. 

(3) The solution to the forced vibration problem is ex- 
pressed as the collection of all possible combinations of the 
plate eigenfunctions and is introduced into the appropriate 
nonhomogeneous differential equation governing the plate motion 
which contains a term representing the surface force acting on 
the plate [the surface force is the sum of the applied force 
and the fluid pressure (in terms of velocity potential expressed 
by the means of spheroidal wave functions)]. The resulting 
equations are then reduced to algebraic equations in the two 
Sets of arbitrary coefficients (plate motion and fluid motion) 
by using the orthogonality property of the "in vacuo" eigen- 
functions of the plate. 

(4) Because the velocities of the plate and of the fluid 
must equal each other at the interface of the plate on fluid, 


the appropriate derivatives of the plate displacement and fluid 
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velocity potential are equated.  Orthogonality of the spheroidal 
wave functions then yields another set of algebraic equations. 
(5) Combining (3) and (4) and solving the result gives 
equations which express the pressure at every point in terms of 
spheroidal wave functions, thus leading to a knowledge of pres- 


sure at the point of interest at a given frequency. 


25 SOLUTION 

The following steps are from the work of Alper [16]. As 
seen in Figure l a circular plate is clamped in an infinite rigid 
baffle, represented by the x, y plane. The upper part is filled 
with fluid and there is vacuum in the lower part. A spatially 
arbitrary, harmonic, transverse force is applied to the "in- 
vacuo" side of the plate. The objective is to calculate the far 
field pressure in the acoustic fluid produced by the forced motion 


Beetne plate. The equation of motion in the fluid is 


ENE ІШЕ) 

where 

£12 ka, frequency parameter 

T-wt, dimensionless time. 

@ = dimensionless velocity potential 
р = 2 90 dimensionless pressure (2) 

AT 

V = EV dimensionless velocity (3) 


EGIT look for harmonic solutions of the form 
^ à 
Ø = D е ic (4) 


where 
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^ 
D = magnitude of velocity potential and is a function of 


space coordinates only. 
Mem@ec, putting Eq. (4) into Eq. (2) gives 

р = ια φ εἰ - pet © (5) 

Again, ind is a function of space coordinates only. 

Pease, putting Eq. (4) into Eq. (1) gives 
ay? (Beir) 0242 d eit - ο 

(ay? +02) fH = ο (6) 

Eq. (6) is the dimensionless Helmholtz equation. Its solutions 
must satisfy the boundary conditions in the fluid. Except for 
the conditions at the plate, these are: 

(a) The fluid velocity normal to the rigid baffle must be 
zero. 

(b) Only outgoing waves are accepted. 

Oblate spheroidal coordinates, which are shown in Figure 2 
[Van Buren, King, Baier, and Hanish, Ref. 17], are the most 
Suitable for these boundary conditions. 

When the confocal ellipses and hyperbolas having an inter- 
focal distance of plate diameter 2a are rotated about the minor 
axis (z axis), flattened ellipsoids and hyperboloids of one 
Sheet are obtained. Under the transformation to dimensionless 
coordinates 7 and br 


a 1(1-97) (%7+1)1172 соз Ө 


к = 
у = а [(1-12) (32+1)11/2 βίπθ (7) 
2 = а i. 


The $ = constant surface is an ellipsoid of rgolution. The 


cross section perpendicular to the x or y axis is an ellipse, 
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ue that perpendicular to the z axis is a circle. In the xy 
plane the eccentricity of the ellipse TRE the semimajor 
axis is $ . Hence the bigger 9 58; the larger the ellipse, 

and the more nearly circular it becomes. The surface? = 0 is a 
Secular disk of radius a in the xy plane. The surface %= 0 

is a sphere of radius ©. The 7) = constant surface is a hyper- 
boloid of revolution of one sheet with anasymptotic cone whose 
generating line passes through the origin and is inclined at the 
angle 7)- cos Y, to the z axis. The cross section perpendicular to 
the z axis is an ellipse, sections perpendicular to the x or y 
ENS are hyrerbolas. The SUR e 0 is the xy plane except 
u he circular disk 3= 0. The surface [η] = l is the z axis. 
ime intersection of Y and ? constant surfaces is a circle. The 
angle 6, measured from the x axis towards the y axis determines 


EB snt on this circle. 


Domains are: 
0647161 , --ς 5 6 ο О <О&?т 


Eq. (23) becomes, in oblate spheroidal coordinates, 
Lee ^ 
ΜΜ νο d Ero uan PME € 8 
(8) 


The general solution of this by separation of variables gives 


A 


Q = 5 Cia ὦ) 5 ia?) OMC) 


(9) 
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where OM С) _) Os md and 


sin MEG 
S n (i, YD and Bes (-10, 13) represent angular and radial 


functions respectively. Putting the general solution Eq. (9) 
into the Helmholtz equation Eq. (8) leads to two ordinary 


differential equations. 





8 d ) m 
feng S Co E 


(10) 


8 а 15) 6, ал) _ A lm 


ІЗ К С52,797-0 (11) 


= 


where A is the separation constant and m = 0,1,2,3.... 


Erter Mathematical calculations 


2 Oo (4) 
Darr = Paya ж (02,7) Kena (1013) 8), 6) (12) 


(n-m) are even values 


where A = constant, 
mn 
оо 
EE (-i2, = > ΛΚ...) 
mn η q m+q 7 
q=0 
MS the oblate angular function of the first kind, and 
m 
"mq (m) 


is the associated Legendre function of order m and degree (m9). 


Eom [ref. l6, pp. 17-23]. 
ПО, —:0, 13) = ain, iS) -i Қы) Cia, 13). 


mn 
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In the development of Eq. (12) some of the possible solutions 
are rejected. The coefficient of the second independent solu- 
Buon of Eg (4), which is σα 10, η) ШЕ τε ο 
S use at T = 1 this function goes infinity. The coef- 


Bent Of the other independent solution combination of Eq. (11), 


which is RC τω, i3) , is also chosen to be zero, because 
only outgoing waves are considered. gt & SE (-ilL, 15) can 


be calculated in power expansion given by Flammer [ref. 18]. 


Also by expanding Eq. (3) in oblate spheroidal coordinates yields 


= | 
V 2 cU = ET 20 5 $ 
A A AAA 9 
(13) 
E. — ЕЕ | 
where е, 7 е, апа е are the unit vectors normal to the 


4, and © surfaces, respectively. The velocity component normal 


to the Y surface Ta 


2 JET ъф 
Vy = 45-5 2X (14) 


The baffle plane corresponds to Y= ο The rigid battle 


condition requires that 


_ 32 (T5) 


V 
E = 





1=0 


Шапо Ба. (12 for Eq. (15) yields 





E (1) | 
dp 6 


(16) 


ы 
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After using the power expansion for San (20,9), odd 
values of (n-m) simply could not satisfy Eq. (16). Therefore 
odd values of (n-m) are rejected and hence it is found that 
Eq. (12) exists only for even values of (n-m). 


Summation of all acceptable solutions yields 


oo Qo 

A . 

$ = > Ann San Ciria 7) O (ө. (17) 
m=0 


n=m,m+2,m+4 


In the case of axisymmetric motion, m = 0 and then the 


far field solution becomes 


oQ 


^ iF Y 
g- = 22—— ЕН a s (-12,7) (18) 
of E On on 

E 


С. MINDLIN-TIMOSHENKO PLATE ANALYSIS [REFS 6 and 16]: 
The dimensionless equations of motion in polar coordinates 


for simple harmonic excitation are: 





(19a) 
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L 
n Co dE | | ο(αθ 
ΣΕ S| or = t alu, )+p, (u, de 
uoc) 
and Q- hp Wa 4 (1) 2s 
gu a LL o9 


uJ(w9)- dimensionless transverse displacement 


СЗ ἵ-ν 


| p 
БЕП CO 
Ve" au? БОЛ προ 
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_h 
v= Yavin 
2° 
Е empirical shear coefficient for plate. 


Y Bana (ρα, θ) an тс ее (Сс косессоп Өс спе nor- 
mal to the plate's middle surface a result of bending in the u 


and O direction respectively, 


ү 


u = — 
a 
πω Ea spatially distributed harmonic external load 


A 
P. Qu, 8)- pelin πας εἴτε pressure to Q (u,6) 


ϕ, - 9 


ος πο ο -τοη) solutions to Eq. (19) are 


A A ^ 
Qu r | oH 
(1) =, κο.) EU e 


^ 
Bor Q 


^ 


Уш 
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^ EN A ^ 
ος cr oe — T SS (20) 


satisfy 





 пЕезлопс to Eq. (21) under the condition that they will 








be finite at u = 0 are: 
u, - 2, επ, αὶ ON 
M=O ie 
^ Р E Y 
Ws, = y ρα IN CN Om <L > yo с^ 
m=0 | Cu 
ü = al ea Im CAD On (22) 
and 
E-. 2 БА J, Ç=, u) On 
м = О a ae 
A о a ας τς 
ο), Ξ > В, С, (11 ш) tA ών 2 
m=O Es? 
H - > 8. Ls) @, (23) 
m=0 
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Since J, (tx) = L, (<) only Eq. (22) can be used without 


losing generality. 


Boundary conditions for a clamped plate are 


ur (1)- 0 
UNDE 0 
Y, (1,9)= 0 (24) 


Eonbining EO. (20) and (22) yields 


gj "EXE pes OS Ja а ра, υ) 8. 8 (6) 


WOSA, ze RER +) u e^ EU | ja 6,9) 
4,6) b i e +, (κιν) En) G) (e) τ 


a 
where ©,- 716 


 Узпа the boundary conditions, Eqs. (24) to (25) and the re- 


quirement that solutions will be non-trivial the characteristic 


* (e; 6; ) Js Cot ) Ja, (9 J C) 


equations are obtained 


Неее) ааган) 


= 0 (26) 

uzi 
where m = 0,1,2, ... and derivatives are w.r.t. u. The dimension- 
less frequencies SL. "EINE OU kes IE M --υ-, τος which 


mj 7 e © è / 


Eq. (26) is satisfied are the eigenfrequencies to which there 
are corresponding values lmj'  2mj' 3mj^ lmj and mi The 
corresponding eigenfunctions are 
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c E Was (u) 8), (е) 


DM Jon (ooi) a)! 
= u C τ. Cu . “, (6) а 
L NM.) 9 Е 


E . Y. к u) G), (6) 
-(6, m 1:7 m Sheng) ы m+ Bo Lep τ Cans αν 





Тала) о б) M 
3 ( S. Е M τ) WA IRG mJ ) E Cr u) (u) (9) (Ὁ) 
m Im amj) = mi n Kim) n 


^ d ζω 
nd Ye PORO- eg EUG {ρα ies Eee s ΜΟΝ 


= (Sins чыте уте d @ (9) (e) 


EL 1,2, ... ; E DE s 


Mm the case of axisymmetric motion, m = 0, Eq. (26) takes 


the form (eliminating the derivatives) 
ad La rt eo PEE" 


and the eigenfunctions 


2] 








ο ο. de Cj) 
Ww J, Cty ) 7 hr) E (Ku u) 


5 E INCON ШЕ | 
J | Ph i) 





= TE - J (eu) + (5; τ΄ 
OB ST ϱ) 3 J 
N 
25 ο (29) 
91 
А 
because @ (9) = 0. Thus, there is no coupling between М/ 
^ ^ 
and in axlsymmeťtric motion. 
W, ana t 
There can be rotations ur 0 that are axisymmetric, but 
N 


Я y Е ^ 
they will be independent of w and Y 
The characteristic equation for independent axisymmetric 


A 
Ше motion is obtained from Eq. (28) by setting the factor 


a) = 0. 30 
x, J, CR) (30) 

When the roots of Eq. (30) are ne into e 2 SH NEP 

values are obtained for axisymmetric Y. motion. The ap- 
lied force is normal to the plate here and these modes are of 
no interest. This can be seen easily after calculating a value 

= 2- 
2 28 nie 
EU ° =i! + O > ) 
| п (оу | ή 
h 

For a thick steel plate, for example, va ΠΡ»; and 
for a thin plate fl would be even greater. 

uumally, Eg. (28) for this problem takes the form 


CS - ne, J (5) 94040 - CS; - 08,1 (0602 


| 4) - 0 (31) 


150% 
By using Eq. (2la) the roots f ӨШ μπιτς equation сап be ob- 


Ἑ ο οἱ and Ξ 0.2 


ше 
pio 


tained.  Symmetric eigenfrequencies for 
Ere given in Table I of Ref. 16. 
2: de 
2224018 | 
«5/0 a Т fv 0 
when у.) (32) 
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Then calling 


e^ t - ос? when ο pHSECHSracterrstic Ес. (31) for 


axisymmetric motion becomes 


T-1)%, TB) J, coh) + (6-1) 7,06) (ËB) = o (33) 


and corresponding eigenfunctions are 


E INC 
We = onu - E" L (Bu) 


^ 
HR Sn = P qe ο Je cos) BT L(g ae 
L. (Bp 


meee COUPLING OF ACOUSTIC FLUID AND PLATE 

The following developments are also taken from Alper and 
Magrab [6]. On the 3 = 0 surface, which represents the piate, 
boundary conditions are: The normal fluid velocity and the 
acoustic pressure on the fluid are equal to the transverse 
plate velocity and the pressure on the boundary of the plate, 
respectively. 

sing Eq. (30), 

EIA E os 
no Aya 9A 1-0 4 25 2 


á 
| 
I 


ον 
Q 
C > 


(36) 


The external force (applied force per unit area plus 


acoustic pressure at the plate boundary of fluid) is given as 


WFT) 


(37) 
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Bence, if Eq. (37) is substituted into Eq. (19c), then, use of 


F (19) and (36) couples the motion of the plate and the 
A ^ 

EEuwstic fluid, so that solutions for Ф and UJ can be found 
т 


under the specified force E a И NDS | - 0 апа 
^ 
ер the eigenfunctions ({АГ,, ‘Ve satisfying 
Bembined Eq. (19) and Eq. (37) could be pr 
Let the solution for the plate motion be expressed as expan- 


mens im its eigenfunctions: 


" ΓΖ». 


MI E 
^ = D ^ 
E rus. 
A NEN AT i. 
^ оо > 
>. Жы To nt T" 
A= O =l 


Putting Eq. (38) into the linear and homogenous combination 


mera. (19) and (37) there is obtained: 


(=) > ж: D eg αφ 


m=0 = 0 
Я > 95 7 ^ 
3] E. = DR AU Di) az: - 0 
MIO Jal 
BS Sp 78 О 
3 > 2. My Car a Las) us E id Е 


The orthogonality condition for a clamped circular Mindlin- 


Timoshenko ο ος is defined as 


E BE ami Y 0 99 T м0 


= 5а 5, (1- Som) T Nj (40) 
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is the Kronecker delta and 


where δις 


\ 
2 d 2 
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Using Eq. (27a) and performing appropriate integrations 
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Lo h 2. 
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For the axisymmetric case, m= Q. 
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‘Multiplying nefasto Ec. (39) by v ᾧ, poa. 
the second by p Pom; udu Je nonsese EE ру 
Wmi u du 49. adding the three equations and integrating over 
0 € O < , OS u$ and using the orthogonality condition 
Eq. (40) leads to 
air | 


2. А: ^ 
8 5,0 ELM (α--ᾳ Ὁ) A] | ES μάμάΘ (42) 


ο ου ο. 5 1,2 
^ 
By expressing ЕК7 іп а Fourier series £(u,@ ) = 


2, FG) €), C8) 


меде 
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[ fue) (9), (6) 46 (43) 





F_ (u) = 


| 
Т 4 TEL ic 


^ 
and using Eq. (12) to express Ф there results 

оо бо ро (4) 
^ 2. 0) | 

Ж I | τ H 

4-7 £6 (9. BE a) ΣΣ B.S. Gant Ga (e) 

m=0 IA NAM ML MTY (44) 
mning Eq. (42), (43) and (44) and performing integration 


with respect to 9 yields 
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BU Эа, а -f W, qu) E E τ; A u du 


(46) 
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MT -{ Fan (a) Won Cu) udu 


Expressions for g . are 
m,nj 


NOR dmarz 2 Kamp 26 


dm nj = 
мат – 5 = MT ) Хм < O s 
where BR = А0807 
Эмате= > om е κ ет 
Je 2 2 о, С. i (m+q)| еман) 2 
C - 1,2 and 


ро = a 
Ga n o o иад 
Π.Ν p 2, q! mme г, 


For the axisymmetric case m = 0. 


Membanina Eq. (12) and (36) and using the first of Eq. (38) 


for w yields 








= 16% 
2,3, W. „w =L5 A Mos 4 - ein a 
el is -- ο, 


(48) 


The functions ` s: _ iM, (n-m) even, are orthogonal 
on the interval © ςηςι ο πο ο οςσομα ву condition is 
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fs (чат) Ж (42,4) = Inn Sn 


(49) 
(n-m) even; m = 0,1,2, 
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E (50) 


where 


(n-m) even. 
 ІГіріуіпс Eq. (48) by Seh ( - iQ, qma, integrating over 


the interval om βάψιμο for An and using the definition 


moe (46) yields: 
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m MU 2n zm. t2, m + 4. 


Econo Eq. (51) into Eq. (45) 
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m = 0, Ер 25 — CESE; 2 





This represents an infinite number of sets of infinitely 


long linear non-homogenous equations in the coefficients Dj: 


mer each value of m, there is a set of equations, t = 1,2, 
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The infinite set of equations, yields solutions for Daj’ 


2, ... Which, when put Еа. (38) gives the displacements of the 


A 


plate and, when put into Eq. (51), gives the coefficients Ann” 
th harmonic 
en 


The coefficients An are then used to determine the m 


of the potential field using Eq. (12) and thereby the m har- 
monic of the pressure using Eq. (5). 

Instead of using an infinite set of equations there is a 
necessity to use truncation with a resulting limitation on 
accuracy (N = number of terms in the truncated solutions). Alper 


ШЕ! cOund that for the far field pressure using N = ll for thin 


E = a h 
E = 0.04) and N= 16 for thick x = 0.2) plates gives solutions 


which are accurate to within 0.2 percent. 

Alper calculated patterns for the first three resonant fre- 
quencies for thin E A CEE СЕ 0.2) steel plates. 
These are presented in Figures 3 and 4. Calculated pressures 
in the far field on the plate axis as a function of dimension- 
less frequency for thin and thick plates are presented in Figures 


5 and 6. Physical properties for steel plates with water as the 
C 


E 
fluid medium are: k = 0.86, y = оа 
τ Es 


= 0.42, Y = 0.132 
апа P. = kr 

As a result it was decided that the best choice for building 
a transducer to meet the desired specifications would be a thick 
mate. A sketch illustrating the form of plate displacement for 
Mie τες (4122 1.9) and second (LL = 6.9) symmetric modes is 
presented in Figure 7. 

The computer programs used for the theoretical calculations, 


which were designed for a CDC 6400 computer, were obtained from 
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Dr. Alper. These have been modified for the CDC 6500 and are 
listed in Appendix A. A shortage of time prohibited actual use 


of these modified programs for design work. 
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μπι ΕΕ ΕΕ ΗΡΙ EROCEDURES 


An experimental test of the calculations of Alper and 
Magrab was conducted. А small, clamped edge plate transducer 
was constructed with dimensions which corresponded to the pub- 
lished calculated results for the thick plate, as shown in Fig- 
ure 8. The material used is steel. The ratio of plate thick- 
ness to radius, h/a, was chosen to be 0.2. The thickness was 
chosen so that the second, in-fluid resonance mode would be about 
75kHz. This corresponds to a value of the parameter ιν 6 9. 
Exterior dimensions were determined by the positioning and 
mounting hole on the cylindrically-shaped torpedo extension sec- 
tion, which served as the acoustic baffle for the transducer. 

The height of the outer cylindrical section at the rear was 
 ЕЕеп to be about 3/4 wave length in order to create a high 
acoustic impedance for the plate edge s the neighborhood of 
Hu kHz. 

The driver for this transducer consisted of a piezoelectric 
ceramic longitudinal vibrator, similar to the sketch shown in 
Figure 9. The ceramic, a lead~-zirconate titanate type, con- 
Sisted of two axially polarized circular discs with a hole in the 
center. Masses on the ends and a clamping screw were constructed 
and their dimensions were selected to make the half-wavelength 
(free-free) longitudinal resonance to occur near 75kHz. A 
Silicone grease was used between joints to ensure good acoustic 


ЕСЕР упс. 
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In order to provide a more nearly point source of excitation, 
the area of contact of the driver with the plate was reduced 
by the use of a small nut on the clamping screw between the 
driver and the plate. Later, measurements were made with the 
full face of the vibrator clamped to the plate. 

Resonance frequencies of the completed transducer assembly 
were determined using measurements of the electrical admittance 
of the piezoelectric elements. A Dranetz Impedance-Admittance 
meter was used for this purpose. Typical graphical results of 
conductance versus frequency are shown in Figures 12 through 14. 

The construction of the longitudinal vibrator used in the 
acoustic measurements prohibited the obtaining of a satisfactory 
isolation of the electrodes from electrical ground of the system, 
needed for making electrical admittance measurements with the 
unit in water. Therefore, resonance frequencies were not mea- 
sured in water. However, measurements were conducted of reson- 
ances of the driver with the transducer in air and in water 
using a different, isolated driver. This driver was a radially 
Eu rrzed Cylinder of barium titanate, 2.54 cm long by 1.27 cm 
outside diameter. It was clamped to the plate using a screw 
and an end mass. Although the free-free half wave length 
resonance frequency for this element differed by several kilo- 
hertz from that of the other driver, the observed resonances 
of the combined systems differed only slightly. The change in 
resonance frequency caused by the immersion in water was 
within the uncertainty of frequency setting. The effects of 
radiation damping on the sharpness resonance was noted. 
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The acoustical measurements consisted of measuring the far 
field acoustic pressure as a function of angle. The arrange- 
ment of apparatus is shown in schematic form in Figures 10-11. 
A pulse method was used. Although the walls and floor of the 
water tanks are essentially anechoic, the reflection from the 
water-air surface can be significant, necessitating the use of 
a box-car integrator. The gate of the box-car integrator was 
adjusted so that only the sound wave traveling by a direct 
path from the source to the hydrophone was measured. 

BE Eevlindrical baffle consists of an aluminum alloy cast- 
ing about 26.25 cm long and 32.4 cm outside diameter. Plastic 
end caps were used. The wall thickness is not uniform but is 
eur 2.2 cm. A sketch is shown in Figure 12. 

The transducer is bolted into a hole which is located as 
shown. Lead weights were added inside the cylinder to make 
the assembly negatively buoyant. 

The rotation system could be arranged to measure the radia- 
tion pattern about either the longitudinal axis of the cylinder 
or the perpendicular axis. Measurements about both axes were 
made. The mechanical drive used to rotate the transducer in- 
cludes a potentiometer on the gear train which produces a D.C. 
 ГГасе proportional to angle of rotation. This voltage is 
КОО еа to the x axis of x-y plotter. The calibrated LC-10 
hydrophone used as a receiver and the transducer were 
Separated by a distance of 1 m with the hydrophone and the 


EEncer of the transducer-at a depth of 1.92 . The 
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output of the hydrophone was amplified, rectified by the envelope 
detector and fed to the box-car integrator. In order to elimi- 
nate the effects of surface reflections, and minimize the ef- 
fects of standing waves in the tank a pulsed signal was used. The 
box-car integrator, gated to the direct path pulse produces a 
D.C. voltage that is proportional to the average amplitude of 

the received pulse envelope. This D.C. voltage is fed to the 

NE E:buic converter which generates the ordinate of the beam 
pattern. The attenuator is used to calibrate the y axis to 
provide a decibel scale. It is needed because of the non- 
linearity of the envelope detector. 

Fach radiation pattern was calibrated by adjusting the x-y 
recorder and driving voltage to obtain at least half scale 
reading on the integrator, and then by the means of attenuator, 
decibel increments on the y scale were marked. 

Results of radiation patterns about the perpendicular axis 
and the cylindrical axis are shown in Figures 20 and 21 and 
Figures 22 and 23, respectively. These were made with a small 
area of coupling between the driver and the plate. Then, with 
the small brass nut removed from the coupling between driver 
and plate the area of contact was increased to that of the 
disk at the end of the driver. Rotation was about the vertical 
axis. Results are shown in Figures 24 through 26. 

The box-car integrator gate width was suitably maintained 
to ensure that the integration would be done at steady state 
region of received signal. 

After some experimental work was done it was noticed that 
bubbles produced in the water greatly affect the results and 
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reduce loading effects of the fluid. Hence to prevent this 
effect a long brush was used for both hydrophone and transducer. 
Before each pattern was recorded, the source and receiver were 
brushed off in order that no effects due to bubbles would occur. 

A second series of measurements using the composite driver 
with small coupling area was conducted. The objective was two- 
fold. The first was to identify the modes of vibration associ- 
ated with the several resonance frequencies noted previously 
and, particularly, to determine which were the first and second 
modes. The second objective was to recheck the results of earli- 
er measurements. 

For the driver alone (free-free) conductance plotted as a 
wC Elon of frequency is shown in Fig. 13. Two resonances are 
noted. It seemed clear from observing the sensitivity to clamping 
provided by grasping it with one's fingers that the uppermost 
resonance involved flexural plate resonance of the lower and 
mass. The lower resonance appeared to re the usual free-free 
mode of vibration. 

When this driver is mounted on the transducer, a number 
of resonances can be observed. A plot of driver conductance 
as a function of frequency is presented in Figures 14 and 15, 
for the transducer alone and for the case where the transducer 
mounted in the cylindrical baffle. 

A small capacitive displacement probe was used to measure 
the relative displacement amplitude at various positions on 
the outer face of the transducer. The probe consisted essen- 
tially of an 8 mm. diameter flat electrode which could be 
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positioned very close to the plate (about 0.1 mm). This elec- 
trode and the plate itself formed an electrical capacitor. The 
τισ] circuit provided polarization potential for this 
capacitor and measured the amplitude of the alternating poten- 
tial caused by changing capacitance when the plate vibrated. 
Although the probe measures the average effect of change in 
capacitance over a distance of 6 or 8 mm, it was sufficiently 
sensitive to get a reasonable estimate of the relative amplitude 
Of plate motion as a function of distance £rom the center of the 
disk. A schematic representation of the displacement measure- 
ment apparatus is shown in Figure 16. 

A piot of the relative amplitudes is presented in Figures 
17 апа 18 for each of the transducer resonances noted previously. 
The motion for the 01 and 02 modes was re-confirmed after the 


transducer was clamped in the cylindrical baffle. 
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IO BC US LON OF. RESULTS 


ENEENUSBPECTS OF AXIS OF ROTATION ON PATTERN 

Radiation patterns were obtained for the same transducer 
but with different axis of rotation, one parallel to the cylin- 
der axis and the other perpendicular to the cylinder axis. The 
results for two such measurements are presented in Figures 20- 
Mana 22-23. 

It is concluded that although some differences between 
experimental patterns exist, these do not obscure the basic 
features of the radiation pattern to be used for comparison 
between theory and experiment. For purposes of convenience, 
all remaining patterns were made by rotating the baffle about 


its perpendicular axis. 


Peer PECTS OF VARYING AREA OF CONTACT BETWEEN DRIVER AND PLATE 
Some radiation pattern measurements were made in which the 
size of the contact area between the longitudinal vikrates was 
increased from about 5 mm diameter to the diameter of the driver 
end mass, 17 mm. Since the diameter of the driver end mass is 
memoalable to the dimension of the first nodal circle for the 02 
mode of the plate, it should be expected that the details of 
the displacement of the plate at the higher frequencies would 
not be a good approximation to that assumed for the point drive 
case. The radiation patterns made at high frequencies using 
this larger area of coupling were very complicated and difficult 
to interpret and they are not included here. The pattern at 
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low frequencies made with the larger area of coupling were not 
significantly different from those using the small coupling 


area. These are presented in Figures 24 through 26. 


Seu DENT IFICATION OF MODES OF VIBRATION OF THE PLATE 

The existence of a number of resonances in the neighbor- 
hood of frequencies calculated for the 01 апа 02 modes of the 
transducer plate appear to be due to the frequency dependent 
impedance of the various parts attached to the plate. The 
driver element is one such attachment. However, the boundary 
conditions at the plate edge are probably more important. The 
massive cylinder at the rear of the plate was intended to 
provide a high mechanical impedance at the plate edge to simu- 
late the clamped edge condition. An outer flange of thickness 
comparable to (but larger than) that of the plate itself was 
needed to secure the transducer in the baffle. There were 
several modes of vibration which involved significant vibration 
of the outer flange. Some of these modes were damped where 
the transducer was clamped into the baffle. It is concluded 
that the resonances at 31 and 63.9 kHz were associated with 
the clamped edge plate modes having mode number 01 and 02, 


respectively. 


D. RADIATION PATTERNS 
The comparison between the radiation patterns calculated 
by Alper [16] and those measured in these experiments must take 
into account several factors which can give rise to differences. 
First, the assumed clamped edge boundary condition is only 
approximated in the tests. 
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Second, the baffle used is cylindrical in shape and is 
therefore neither infinite nor plane, so the assumption of the 
theory of an infinite plane baffle is not well approximated 
even though the diameter of the cylinder is about 7.5 times 
the diameter of the plate. 

Mard, the fact that the baffle is not perfectly rigid per- 
mits some coupling of flexural wave energy from the plate into 
the baffle and the possibility of additional radiation of sound 
into the water from the baffle itself. 

Fourth, the area of contact between the driver and the center 
of the plate was made as small as possible but it was not a 
point driving force, as assumed in the theory. 

The comparison between theoretical and experimental radia- 
tion patterns should best be made by inspecting Figures 3 and 4 
and 27 and 29. The theoretical patterns due to Aem (shown in 
Figure 19) have been redrawn in the rectangular log-linear plot 


format used for the experimental data. 


The experimental pattern which should correspond most closely 
to the Ol mode is Figure 27 at a frequency of about 3l kHz. 

Those which should correspond most closely to the 02 mode 
are shown in Figures 20-22-29 at a frequency of 63 kHz. 

For the 3l kHz pattern it is seen that there are some small 
Escutlations of about 3 dB near the center of a single major 
lobe. The level of this lobe at the 90? angle is down about 
l8 dB below the axial value. This contrasts with the theoretical 


90° value of about -2 dB. The reasons for this disagreement 
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are not known, but it is believed that an important contribu- 
tor is the finite size and curvature of the baffle. 

For the 63 kHz patterns, there are some features which 
correspond well to the theoretical pattern. The angular 
position of the first minimum is at about 23 degrees, very 
close to the theoretical value of 22.5 degrees. The level of 
the measured maximum of the size lobe is about the same as that 
of the major lobe. The theoretical curves indicate a side 
lobe maximum level of -2.5 dB located at 48° and a level at 
90° of -5.6 dB. The experimentally observed level at 90° is 
about -20 dB. Again, this lower level near °0° is believed to 
be due largely to the finite baffle dimension. The additional 
structure on the measured patterns is believed to be due to 
radiation of sound from the baffle. 

meLustrcations of the complexity of radiation patterns at 
other higher frequencies are shown in Figures 21-23-30. The 
frequencies here correspond to resonances in which there are 
additional displacements of the transducer parts particularly 
the outer flange (see Figures 16-17). Very likely the baffle 
itself was also excited into some flexural modes of vibration. 

It is concluded that the major featuresof the radiation 
patterns measured agree with those of the Alper theory and 
that the differences are probably due to boundary conditions of 


the plate which did not meet those assumed in the theory. 


E. RECOMMENDATIONS 
Since the goal for this effort was to develop the Alper- 
Magrab theory as a design tool for the development of a flexural 
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disk transducer that might be suitable for use as a sound 
source on underwater test range tracking, the following comments 
and recommendations seem to be appropriate. 

The adaptation of Alper's program to fit the CDC 6500 Com- 
puter, HAL, is virtually complete. Shortage of time precluded 
completion of this part of the work. It is recommended that the 
work be finished and that it be used as intended. 

It appears that some sort of variation on the radiation 
pattern from the second radially symmetric plate mode does offer 
the potential for controlling the source level near the axis 
and yet maintain a relatively broad beam width for the secondary 
lobes. The patterns observed hereshow a deep but narrow null 
at about 22 degrees off the central axis. It may be possible 
Meee this null could reduced by proper distribution of excitation 
force on the plate. It also is possible that the effect of 
this narrow null on range tracking is acceptable. It is recom- 
mended that the experience of this work be employed to con- 
struct and test on the tracking range a flexural disk trans- 
ducer which operates in this mode. 

The smoothest radiation pattern is given by the plate 
operating in the first symmetric mode. The transducer used in 
these tests was not optimized for this mode. It is likely 
that a transducer designed to operate in this mode at 75 kHz 
would have such a small plate diameter that a problem would 
exist in fitting an adequate sized driver into the available 
Space. If the program development tool becomes workable, cal- 
culations should be made to test the design at 75kHz to deter- 
Mine feasibility. 
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A much more likely application for a flexural disk 
source operating in the first mode appears to be at frequencies 
significantly lower than 75 kHz, such as are used on other test 


ranges. 
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Figure 2. Oblate Spheroidal Coordinate System 
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Figure 9. Driver Elements and Position on the 


Transducer 
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Figure ll. Radiation Pattern Measurement Apparatus 
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RELATIVE PLATE DISPLACEMENT 
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Figure 17. 
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Figure 18. Relative Plate Displacement for 01 and 02 
- mode (transducer mounted in baffle in air) 
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Figure 19. Theoretical Radiation Patterns of 01 and 
02 Mode. 
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Figure 20. Radiation Pattern (small area, perpendicular 
axis) 
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Figure 21. Radiation Pattern (small area, perpendicular 
axis) 
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Figure 23. Radiation pattern (small area, cylindrical axis) 
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Figure 25. Radiation Pattern (large area, perpendicular axis) 
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Figure 26. Radiation Pattern (large area, perpendicular axis) 
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APPENDIX A 


COMPUTER ReEReGraMs FOR FAR FIELD PRESSURE 


AND ASSOCIATED PARAMETERS 


Alper's original programs were not changed very much. 
For the main program, the outputs of the Program, Roots, and 
the Program, Caller, (for getting sum)are needed. After 
getting the output of Program, Roots, Program, Caller output 
is obtained by using these results: Finally, the main pro- 
gram will be run by using the obtained parameters and sums. 
For this reason some Formats were changed to fit inputs and 
outputs each other. The only problem which still was not solved 


a8 that of Subroutine Bessel. 
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